Locating the source of projectile fluid droplets 
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The ill-posed projectile problem of finding the source height from spattered droplets of viscous fluid is a 
longstanding obstacle to accident reconstruction and crime scene analysis. It is widely known how to infer 
the impact angle of droplets on a surface from the elongation of their impact profiles. However, the lack of 
velocity information makes finding the height of the origin from the impact position and angle of individual 
drops not possible. From aggregate statistics of the spatter and basic equations of projectile motion, we 
introduce a reciprocal correlation plot that is effective when the polar launch angle is concentrated in a 
narrow range. The vertical coordinate depends on the orientation of the spattered surface, and equals the 
tangent of the impact angle for a level surface. When the horizontal plot coordinate is twice the reciprocal of 
the impact distance, we can infer the source height as the slope of the data points in the reciprocal correlation 
plot. If the distribution of launch angles is not narrow, failure of the method is evident in the lack of linear 
correlation. We perform a number of experimental trials, as well as numerical calculations and show that 
the height estimate is insensitive to aerodynamic drag. Besides its possible relevance for crime investigation, 
reciprocal-plot analysis of spatter may find application to volcanism and other topics and is most immediately 
applicable for undergraduate science and engineering students in the context of crime-scene analysis. 
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I. INTRODUCTION 

The impact of spattered droplets of viscous fluicP^ on 
a horizontal surface results in elongated impact profiles, 
from which it is easy to locate the vertical axis of ori- 
gin from the orientation of elongated impact profiles (see 
Fig. 1). It is accurate for blood-like fluids (and routine 
practice in forensics 3 -* 4 ^ to estimate the angle of impact 
of a viscous droplet from its impact profile as 



sin 



profile width \ 
length / 



(i) 



The subscript "I" denotes impact-related quantities. 
This construction, which assigns to each elliptical pro- 
file (see Fig. 1, inset) the projected shape of the incom- 
ing spherical droplet on the surface, works remarkably 
well. From estimates of 1 from Eq. ([!]) we would like 
to estimate the location of the droplet launch from the 
positions and angles of the impacts. 

However, even when its vertical axis is known, the 
height of a source is not deducible from the location and 
impact angle of individual drops {r l and I in Fig. 1). 
The impact velocity of each drop is missing the infor- 
mation needed to infer the source height, and thus an 
height must involve assumptions. A method widely used 
in forensics 3 is to extrapolate straight lines from each im- 
pact back to the launch axis and seek a minimum in these 
height values, which is consistent only with the equations 
of projectile motion in the limit of high velocity and low 
aerodynamic drag. Otherwise, the strategy systemati- 
cally overestimates the source heiglrP (see Fig. 2). Such 
considerations have led to alternative proposals, such as 
seeking missing velocity information in the detailed struc- 
ture of impact profiles. 4 

In this paper we present a statistical and graphical 




FIG. 1. Profile orientations can be extrapolated (dotted 
lines) to a horizontal point of convergence (o) so that the ver- 
tical axis of origin and the distance ri are known. Even with 
knowledge of the impact angle 0i, we still cannot determine 
the source height (for example, zo or z' ). (Inset) Digitized 
photograph of experimental impact profile fit to an ellipse. 
The impact angle is obtainable from the residue profile width 
and length by attributing the elliptical profile to the projected 
shape of the incoming spherical droplet. 



method of back-estimation consistent with the equations 
of projectile motion and provide a height estimate in 
cases where the droplets are launched within a narrow 
range of the launch angle, using only the impact location 
and inferred impact angle. In the simplest case, linear- 
ity appears in a plot of the tangent of the impact angle 
versus twice the reciprocal distance of impact. 
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FIG. 2. Position and impact-angle dat a (o) for floor spat- 
ter generated as described in Section IV (clapper horizontal; 
i.e. launch angle 0°). Extrapolated points on the axis of con- 
vergence (0) fail to estimate the source position (•) (height 
88.0cm) even as a lower bound (a high- velocity assumption). 



II. RECIPROCAL PLOT FOR FLOOR SPATTER 

The kinematic constant-acceleration equations^! for a 
projectile droplet in cylindrical coordinates z and r (ne- 
glecting aerodynamic drag) are 



z(t) = z + \ (y z + v z $)t, 



(2) 



where zo is the actual launch height. Along the trajectory 
the horizontal position is r(t) = v r t with v r the radial 
velocity. The vertical velocity is v z — v r tan with up- 
ward from the horizontal. If we substitute v z = v r tan 
and v Zi o = v r tan#o into Eq. (J2|, we find 



tan^ = (zq — 2q) htan^o- 



(3) 



We have defined l downward (that is, as —0) so that l 
is positive at impact as in Fig. 1. For launch and impact 
at equal height (zq = z T ) Eq. ([3| correctly predicts that 
#o =0i- Consider impacts on a surface at z l = (a floor), 
where we have located the vertical axis of the origin (see 
Fig. 1) and wish to find the launch height zq. In Fig. 3(a) 
we plot the quantity for actual spatter data, 



t i (floor) = tan#! 



(4) 



versus the quantity 2/r l . According to Eq. Q our data 
will satisfy 



tan 0q . 



(5) 



The slope, that is, the coefficient of 2/r l5 gives the launch 
height zq. Figure 3(a) shows a trial using a standard 
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FIG. 3. (a) The floor spatter data of Fig. 2 in a reciprocal 
correlation plot, using T\ — tan^i from Eq. Q. The slope 
yields an estimated launch height of 95.0 =b 2.2 cm and an es- 
timated launch angle of —5.8 =b 2.4°. Pearson's coefficient of 
linear correlation is r = 0.996. The actual source height is 
88.0 cm. (b) Trial with launch height 85.0 cm with the clap- 
per sideways for a maximally broad launch angle distribution. 
In this case the reciprocal plot yields no height estimate, as 
evidenced by lack of linear correlation (r = 0.143). 



"clapper" mechanism (see Sec. IV) in which a reciprocal 
correlation plot yields the source height as 95.0 ± 2.2 cm 
compared to the actual height of 88.0 cm. The remaining 
discrepancy may indicate residual systematic error in our 
procedure, such as a surface effect on profile shape or 
some other factor. 

To further illustrate the predictions of Eq. (J5|, 
Fig. 4(a) shows how a varying source height above a level 
surface leads to different slopes in a reciprocal correlation 
plot. Figure 4(b) shows how varying the launch angle of 
spatter leads to different vertical intercepts. 

Fig. 5 applies a reciprocal plot to spatter data col- 
lected from walls and ceiling, as well as from the floor. 
Extensions of Eq. Q necessary for this are discussed in 
Sec. El 
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FIG. 4. (a) Reciprocal correlation plots for floor impacts with 
horizontal launch: (•), launch height zo — 88 cm (the same 
data as in Fig. 2(b)); (o), launch height zo = 164 cm (o). 
The clapper orientation is horizontal for both trials. Note 
the differing slopes, due to their differing heights, but similar 
intercepts. (Upper curve: slope 164 ± 6 cm, intercept Oo = 
2.7 ± 3.7°, r = 0.996. Lower curve as in Fig. 3(a).) (b) 
Plots for floor impacts, with the clapper at Oo — 45° (•) and 
Oo = 0° (o), same data as in (a)). The heights of 165cm and 
164 cm, respectively, are similar. Note the similar slopes, but 
distinct intercepts. (Upper curve: slope 168 ±4 cm, intercept 
0o = 45.8 =b 1.2°, and r = 0.995. Lower curve as in part (a).) 



III. AERODYNAMIC DRAG 

In reciprocal plots such as Figs. 3(a) and 4, points 
near the plot origin represent distant impacts, with high 
launch speed. Here an error due to aerodynamic drag 
might be anticipated. Our trials suggest no obvious prob- 
lem, and therefore we turn to numerical calculations. The 
drag force on a droplet with speed v i d 1 | 2 | T | 



Fd = ^PairAv C d , 



(6) 



where A is the cross-sectional area, p a i r is the density of 
air, and Cd is an empirical coefficient. Separate calcu- 
lations show that for our spatter trials Cd ~ 0.5. We 
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FIG. 5. Reciprocal correlation plot of the combined data from 
the floor (o), walls (•), and ceiling (<0) with t\ generalized 
as described in Sec. |IV[ For the combined data the slope is 
zo = 169 =b 7cm and the intercept is 35 =b 7° (dotted line), 
compared to the actual values of 165 cm and 45°, with r = 
0.973. The separate slopes and intecepts of the floor data 
(height 159±21cm and angle 48±13°), wall data (167±18cm 
and 34 ±14°), and ceiling data (185 ±14 cm and 14 ±9°) are 
not statistically distinguishable. 



define 



and 



/ 4 pug gd 

3 Pair Cd 



9 



(7) 



(8) 



where g is the acceleration of gravity, v T is the terminal 
velocity of a drop of diameter d and density pn q , and z T 
is a "terminal height" over which the terminal velocity is 
nearly reached in a vertical fall (y = 0.93 v T at z T ). We 
write equations for the drop trajectories which incorpo- 
rate Eq. 



dvg 

dt 

dv. 



dt 



(9a) 
(9b) 



where v = (v% 
variables 



Vy) 1 / 2 and we use the dimensionless 



u x,y 
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(10) 



We numerically solved 9 Eq. (9| and constructed recipro- 
cal plots in the presence of drag, at three launch angles 
(see Fig. 4). These dimensionless plots, rescaled to the 
appropriate v T and z T , apply to any droplet diameter or 
launch height. In all cases distortion is greatest on the 
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FIG. 6. Reciprocal plots for trajectories in the presence of 
aerodynamic drag (•) calculated using Eq. ([9]) at launch an- 
gles 0°, 22.5°, and 45°. The source height is zo = 0.1z T and 
f i = ri/z T , where z T is the characteristic fall height leading to 
the terminal velocity. The horizontal axis is rescaled by 10 -3 
for rough correspondence to cm -1 in our experiments. In the 
linear region the slopes are little changed from the straight 
lines (dotted) which would occur in the absence of drag. 



left side of the plot, as expected, where the range and 
launch speed are highest. 

A notable feature of Fig. 6 is that, for moderate height 
and launch speed where the plot is linear, aerodynamic 
drag mainly affects intercepts, and thus the inferred 
launch angles, but not the slopes. The height of ori- 
gin estimates are only weakly affected. This dependence 
indicates a certain robustness of the method in the pres- 
ence of drag, as long as we obtain linear plots, which 
might help explain the success of our plots in Figs. 3(a), 
4, and 5. 



IV. METHODS AND STATISTICS 

In the spatter trials whose results are plotted in 
Figs. 3-5, a viscous fluid was spattered on floors, ceil- 
ings, and walls with a clapper device modeled after those 
used in forensics training. 3 Two wooden boards, joined 
at the rear by a spring- loaded door hinge, slammed shut 
at the front, an impact area that we fitted with metal 
plates. A small pouch of fluid to be spattered was taped 
to one plate. The results in Figs. 3(a) and 4(a) were 
obtained with the impact plates horizontal. To vary the 
launch angle we tilted the clapper forward and backward. 
For a broad distribution of launch angles we rotated the 
clapper about the long axis so that the launch area was 
vertically oriented. 

We found that reliable impact angle estimates required 



considerable practice, including self-calibration using ver- 
tical drop impacts on inclined surfaces, during which we 
found that actual blood drops resulted in well-defined el- 
lipses and surprisingly good impact angle estimates using 
Eq. g. To avoid using actual blood in our spatter tri- 
als, we used a heterogeneous fluid whose droplet impact 
profiles closely resemble those of blood droplets. The ma- 
jority of our experiments used approximately 2-3 parts 
Ashanti chicken wing sauce (Bridge Foods) with 1 part 
Ivory dish soap to aid cleanup, and trace amounts of food 
coloring to enrich the color for digitization. In other ex- 
periments, a viscous test fluid of 4 parts corn syrup to 
1 part water with food coloring also led to a reasonably 
successful analysis. 

From a typical spatter trial we chose approximately 
twenty profiles having well-defined elliptical profiles, rep- 
resenting a range of distances. For the plots, we took 
digital photographs of the profiles, and used image anal- 
ysis software 8 to manually fit the ellipses to the profile 
outlines shown in Fig. 1 (inset). We obtained best results 
by matching well-defined portions of the elliptical outline. 
Very close to the source, near-circular profiles can lead 
to large uncertainties both in Eq. ([!]) and in 2/rj when 
r l is very small, 4 and hence some of these points were 
discarded. 

Our estimates of slopes, intercepts, and their errors 
were made using standard regression techniques, 9 where, 
for example, the estimate of the slope (that is, launch 
height) is 



^est — ^0 



2z 



(11) 



with zq equal to the actual launch height and Au = u — 
(u), where (• • • ) denotes an average. A widely known 
parameter of linear correlation 9 is Pearson's r, which for 
our reciprocal plots is 



v /((A^)2)((Ar I ) 2 ) 



(12) 



A complete description of projectile motion requires that 
Eq. Q be combined with the relation 



tan 2 6>! = (r7 2 + 1) tan 2 O + rj 2 



(13) 



which parameterizes lines of constant 0o in the reciprocal 
plot by the variable r] = y/2gz$/vQ, where is the launch 
speed. 

Consider the special case of a product probability dis- 
tribution for the launch angle and velocity P(6q,v) = 
Pi(0o)p2(v), where Pi(#o) is symmetrical about the hor- 
izontal. We can obtain, using Eqs. (|3|, (n), (12]), and 
(13) the simple relation 



2 Zq [for symmetric Pi(0q) 



(14) 



Although Eq. (14) is not general, it suggests that bias 
may occur in height estimates if r 2 is not reasonably 
close to 1. A requirement of strong linear correlation 
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also accords with common sense. Figure 3(b) showed 
an unfavorable trial in which the spread in launch angle 
is maximized by turning the clapper sideways. This trial 
was performed ten times, yielding the r values 0.14, 0.18, 
0.78, 0.86, 0.80, 0.56, 0.83, 0.58, 0.72, and 0.37. These 
values are all smaller than 1, and the plots were uncor- 
rected in appearance. In contrast, in Figs. 3(a) and 4 
r = 0.995 or higher, and for the combined data sets in 
Fig. 5, r = 0.973. 

We can also combine and analyze spatter data from 
floors and other non-horizontal surfaces in a single recip- 
rocal plot. Co-plotting data from non- horizontal surfaces 
(see Fig. 5) requires corresponding redefinitions of the 
variable r Y . We then seek a slope via Eq. ([5| as before. 
Spatter on a ceiling at height z = z Y is incorporated by 
adding a term to Eq. Q and changing the sign of the 
tangent, 



^(ceiling) = z l tan^. 



(15) 



Ceiling data tended to be more robust closer to the 
source. On walls the right-hand side of Eq. ([!]) yields 
the angle W of impact with the plane of the wall. The 
correct vertical angle is then l = sin -1 (cos cos 6 W ) 
with aj the angle from vertical of the wall profile. We 
plot 



t i (walls) 



tan I? 



(16) 



where z l is the elevation of the impact. We typically 
avoided near-horizontal wall impacts. Once the axis of 
origin is located, tan^ in Eq. (16) can be replaced by 
cos^/ tana x , where fa is the angle between the wall and 
the vertical plane of impact, a relation which is conve- 
nient for curved vertical surfaces such as pipes. 



V. DISCUSSION 

We have introduced and illustrated with several ex- 
periments, a plot-based method for locating the spatial 
source of spattered viscous fluid. Th method is effective 
when the spatter is launched within a narrow range of 
polar angles. We showed how a reciprocal plot of the im- 
pact data, together with elementary projectile physics, 
can exhibit linear trends among the data points. From 
the slope of any strong linear correlation that occurs, we 
obtain a robust estimate of the height of the origin that 
would otherwise be unavailable. Broad distributions of 
launch angle cause the method to fail, and in such cases 
we reach a null conclusion (a lack of linear correlation), 
rather than an erroneous estimate of height of origin. 

We also showed that the method appears insensitive 
to aerodynamic drag effects within some velocity regime. 



We also extended the plotting method to some other sim- 
ple geometries. 

The reciprocal plot, which is based on the correct 
equations of projectile motion, may eventually become 
a useful tool for forensic spatter analysis, where informa- 
tion on height of spatter origin pertains, for example, to 
whether a victim was sitting or standing. 

Another possible use of this analysis, which we are pur- 
suing separately, may be its application to geophysics and 
volcanism, especially to phenomena such as volcanic ejec- 
tion and lava fountains.^ 

The most immediate use of reciprocal-plot analysis is 
as an undergraduate activity for science and engineer- 
ing students. In the simplest case students might use 
a wooden block or book to spatter a puddle of viscous 
"blood" (such as corn syrup and food coloring) from an 
elevated surface such as a table or shelf, tabulate the po- 
sition and eccentricities of droplet profiles on paper taped 
to the floor, and analyze this data via reciprocal plots in 
a spreadsheet program to reproduce our Fig. 3(a). Such 
an activity is a novel use of the equations of projectile 
motion within the exciting context of crime-scene inves- 
tigation. 
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